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Abstract 

We consider Toda flows induced on the set of 
orthogonal upper Hessenberg matrices. The 
explicit formulas for the evolution of Schur pa- 
rameters are given. 

1 Introduction 

Any symmetric nonnegative definite Toeplitz 
matrix Tn+i of order n + 1 can be modeled 
as an autocorrelation matrix of a stationary 
signal [12] 

p 

Xm = XI O/ cos(t71w/ + 0t)+ ym, 

1=1 

where 0/ are arbitrary phase shifts and z/m is a 
zero mean white noise process whose variance 
equals the smallest eigenvalue of Tn+i. 

Assume that the eigenvalue Amin is simple, and 
let (r;o, • • • , r)n) be a corresponding eigenvector. 
Then [12] the polynomial 

= Vo + * • • + 

hcis n distinct roots Aj, • • • , A„ on the unit cir- 
cle, and the frequencies of are given by 
{exp(±ici;/}|Lj = {Aj}”_j, where i denotes the 
imaginary unit. One can construct [2] an or- 
thogonal Hessenberg matrix O with character- 
istic polynomial proportional to 'tpn- Moreover, 



the amplitudes q/ can be recovered from the 
first components of the normalized eigenvec- 
tors of 0. One can then use any of several al- 
gorithms designed for unitary and orthogonal 
Hessenberg eigenproblems [9, 6, 1, 10, 4] 
to calculate the frequencies and amplitudes. 
In the present paper we investigate another 
aspect of orthogonal Hessenberg matrices. 
Namely, we consider Toda flows on these ma- 
trices (referred to as Schur flows in [3]) and 
obtain explicit formulas for the evolution of 
the so-called Schur parameters under the Toda 
flow. Since Schur parameters determine or- 
thogonal Hessenberg matrices uniquely, we ac- 
tually obtain an explicit description of the evo- 
lution of a given orthogonal Hessenberg matrix 
under the Toda flow. 

2 Inverse problem for orthog- 
onal Hessenberg matrices 

Let M be the set of positive Borel measures 
on C which have the following properties. For 
any ji £ M the support of /x, which we denote 
by A^, consists of exactly n points which lie 
on the unit circle U and is such that: i) if 
A E A^, then the complex conjugate A is also 
in A^ and /x{A) = /x{A); ii) /x{A} > 0 for any A 
in A^; Hi) fi(C) = 1; iv) —1 ^ We further 
introduce a class OH^ of orthogonal matrices 
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0 =11 oij II such that Oij = 0 if i — j > 1 , 
Ot>i,t > 0 for all i and det 0 = 1. Finally, given 
a vector r = (tq, • • • , G introduce a 

corresponding Toeplitz matrix T{r) =|| t{j ||, 
where Uj = 

Theorem 2.1 Given a positive definite sym- 
metric Toeplitz matrix T{r) with tq = 1 there 
exist exactly one measure p ^ M and exactly 
one 0 G OH^ such that 



f AV/x(A) = r,- =< ei,0*ei >, (2.1) 

Jc 

t = 0, • • • , n — 1 . Here ei , • • • , e,i is the canoni- 
cal basis in and <, > is the standard scalar 
product. Conversely^ for any p £ M and any 
0 G 0^+ the matrices T(r), T(r') are Posi- 
tive definite Toeplitz matrices. Here 

Ti = / =< ei,0‘ei >, 

Jc 

t = 0,- • - ,n — 1. 



Remark 2.2 Theorem 2.1 is more or less 
known to the experts (see e.g. [8], [11]). We 
nevertheless give an independent proof to clar- 
ify relationships between introduced objects. 

Remark 2.3 There is nothing mysterious 
about the number —I which we have excluded 
from the support of each measure in M , This 
simplifies notations a little bit. 



We need the following elementary lemma. 

Lemma 2.4 Let vi, • • • , be an orthonor- 
mal system of vectors in R^, There exists ex- 
actly one orthogonal matrix 0 such that Oe{ = 
t;,-, i = 1, * • • , n — 1 and detO = 1. 



We can now outline a proof of Theorem 2.1. 
Proof: Denote by the vector space of real 
polynomials of degree less or equal n — I, Set 



<A%A^ >= / y-^dp{X), 
J c 



( 2 . 2 ) 



We prove that (2.2) defines a positive definite 
scalar product on P„. Observe that 



Jx'X^dn = J X'-^dfx = jx'X^dn. 



Indeed, f X'X^dfi = _Ea€A^ = 

^A6A^ A‘~-'/r{A}, since A = A“^ for X £ U. 
Further, since /r{A} = /r{A} we have 

A‘-V{A) = ^ A‘-V{A} = / A‘Vd/r. 

A6A,, A6A,, 

Let g = oo + . . . + a„_iA”~^ G P„. We have 

< q,q >= aiUj / X'X^dfi = / \ q dfx > 0. 

t,i=o 

Further, / | g p d/r = 0 if and only if g(A) = 0 
for any A G A^. Since degg < n = card(Aft), 
this is possible only if g = 0. Consider the 
polynomial ^(A) = -t) = bo + ...+ 

h„_iA”~* + A". Since all roots of ^ lie on 
the unit circle we clearly have A”^(l/A) = 
bo^{X),bo = ±1. Further, all coefficients of 
^ are real because = A^. Consider the 
linear operator 0 : Pn -* Pn defined as fol- 
lows: OA’ = X'+\i = 0,---,n-2, OA"-* = 
—bo — 6i A — ... — 6„_i A”~^ . We now prove that 

0 is orthogonal relative to the scalar product 
<, > . We should prove that 

< OA’,A-' >=< A‘,Cr*A^ > 

for any t,j = 0, •••,n — 1. The only non- 
trivial case is i = n — I, j = 0. We have 

< OA""\l >= -bo - 6 iTi - ... - bn-lTn-l, 
where r,- = Jq X'dfi{X). Let 0“^1 = cq + . . . -b 
Cn-iA”~*. Then < A”“^,0~^1 >= cor„_i -f 
CiTVi- 2+- • •+Cn-i. Thus, it is sufficient to prove 
that Ci = —bn-i-i,i = 0, — 1. We clearly 
have 1 = cqOI - 1 -. . .-bCn_iOA"“^ = cqA- 1- . . .-f 
Cji— 2A*^ ^"bCn— i( bo~biX — ...— bn—\X ) Or 

1 = — C„_i6o, Co — C„_i6i = 0, Cl — Cn-lf>2 = 
O,-*-, c„_2 - C„_i6n-1 = 0. This yields 6] = 

— Co/f>0, ^2 = ~Ci/f>0 ,’”» f>n-l — — Cn-2/^0- 

We now use the relation A”^(l/A) = 6o^(A). It 
follows that bn-i = bob{, i = 0, Thus 

bn— i/bo — c,*_i/hoi ^ ” I?***?^' These 
are exactly the required conditions. Thus 
we have constructed an orthogonal operator 
0 such that Jq X'dfi =< 1,0*1 >, i = 
0, 1. Observe that the characteristic 
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polynomial of O coincides with Thus the 
spectrum of O is In particular, detO = 
1 (here we use the assumption that —1 ^ 
A^). Let po = be an orthonor- 

mal basis in Pn obtained by the orthonor- 
malization of the basis 1, A,* • • , It is 

clear that the matrix O of the operator O 
is upper Hessenberg in this basis. Moreover, 
the entries are all nonzero (otherwise, 

span(po) • • • = span(l, • • • , A*”^) is an 

invariant subspace of O which is not true). 
Without loss of generality one can suppose 
that > 0 for all i. Otherwise one can 

take diag((i, - ' • ,(n)Odiag{(i, - ■ ■ ,en). 

Suppose we are given a positive definite 
Toeplitz matrix T(t) and an orthogonal ma- 
trix O € OH+ such that Ti =< ei,0'ei >,i = 
0 , • • • , n — 1 . Then 

T{t) = yTy, (2.3) 

where V is the upper triangular matrix 
[ei, Oci, • • • , with positive entries on 

the main diagonal. But (2.3) is the Cholesky 
decomposition of T(r). Hence it is uniquely 
defined by T(r). In other words, the vectors 
Oei, • • are uniquely defined by T(r). 

Since these vectors form a basis, the vectors 
Oei, • • are uniquely defined by our 

Toeplitz matrix. Thus by Lemma 2.4 the ma- 
trix 0 is uniquely defined by T(r). Given a 
positive definite Toeplitz matrix T{r) we can 
endow with a scalar product <,> and the 
shift operator defined on span(l. A, • • • , A”~^) 
a 5 we did before. Then using Lemma 2.4 we 
can extend this operator to the orthogonal op- 
erator 0 , defined on Pn such that detO = 1 . 
Then the matrix of O in the basis obtained by 
orthonormalization of the basis 1 ,A,---,A”~^ 
belongs to OH^ and r, =< ei, 0 *ei >,i = 
0 ,***,n — 1 . Consider now the rational func- 
tion 

/(z)=<l,(z/- 0 )-'l>. 

As is easily seen 



where all r, > 0 . We then can define the mea- 
sure p £ M by the conditions p{X{} = r, and 
equal to zero otherwise. We immediately see 
that equations (2.1) are satisfied. It remains to 
prove that the measure p is defined uniquely 
by conditions ( 2 . 1 ). Let pi^ £ M^k = 1,2 be 
such that 

/ X^dpi = / A*d/i2, 

Jc Jc 

i = 0 , ***,n — 1 . Then we can construct 
Ok^k = 1,2 such that conditions ( 2 . 1 ) are 
satisfied. But then 0\ = In particular, 
A^i = A^ 2 , i.e.. Pi = p 2 because we have 
for p{X} the following system of Vandermonde 
equations: 

I = 0 ,* •• ,n - 1 . ■ 

Let T{t) be a positive definite n x n Toeplitz 
matrix and <,> be the corresponding scalar 
product on Pn- Let 

Px{X) = <5, A* + > 0,z = 0, •••,n - 1, 

be the basis obtained by the orthonormaliza- 
tion procedure from the basis 1 , A, • • • , A”~^. 
Since pi 

is orthogonal to span(l. A, • • • , A*~^), we have: 
Apt(A) is orthogonal to 5pan(A, • • • , A*). Fur- 
ther, r = Xpi{X)/6i - pi^i/6i^i £ P.+i. Let 

e Pi^i be such that < q,(fi >= q{0) for 
any q 6 P,+i. Since p,- is orthogonal to P,- and 
both r and (pi are orthogonal to AP,, we obtain 

^PiW/^i = Pi+li^)/^i+l + 7t¥’«5 (2-4) 

for some real 7 ,, i = 0, • • • , n — 2. An easy cal- 
culation shows that ipi = ^,A*p,(l/A). Hence 

1 = ^W+i+7M^ (2.5) 

In other words, if we know 7 o,*-*, 7 n- 2 ? we 
can find (5i, • • • , 6n-i- Then using (2.4), one can 
determine pi,* * *,Pn-i and consequently using 
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again (2.4) the corresponding upper Hessen- 
berg orthogonal matrix O. We have by (2.4) 

^ '^/^i(^)? /^t(^) “ 0, • • • , n 2. 

Evaluating (2.4) at 0, we obtain Pt*+i(0) = 
i = 0,***,n-2. Thus =< 

pK*^) ^ 2. 

Further, Oi^i = — 7oPo(0) = —70. Let us set 
~ — l/^M “ 7i — l^t — 1’ 

i = . . .^ n — 1. We obviously have 

erf + = 1 , Oi^i = 

i - i/Q = 1- Further, On^n = 

±yjl — CTn-i- The sign is defined by the condi- 
tion detO = 1. The quantities Ui^Ci are called 
Schur parameters and auxiliary Schur param- 
eters, respectively. As we saw above the Schur 
parameters //,, i = l,**-,n — 1, determine O 
uniquely. 

On the other hand, if we know the entries 
^ = 1 • •• ,n — 1 of the matrix 
O we can determine 7^ by (2.5) up to a sign. In 
other words, the entries (auxiliary Schur 
parameters) determine O almost uniquely. 

3 Explicit formulas for the 
evolution of auxiliary Schur 
parameters under the Toda 
flow 



uniquely defined by its auxiliary Schur param- 
eters (Ti{t) = We now describe ex- 

plicitly how these parameters evolve under the 
Toda flow. 



Theorem 3.1 



£T.(<) 



Aj(t) 






i = l,*”,n — 1, Ao = 1. Here Ai{t) is 
the i—th principal minor of the matrix r{t) = 
exp((O(0) + O(0)'^)0. 



Proof: We 

know [7] that 0{t) = R{t)Q{0)R{t)~^ , where 
exp(O(0)t) = Q{t)R{t), Q{t) is orthogonal, 
and R{t) is an upper triangular matrix with 
positive entries on the main diagonal. We then 
clearly have 



<Ti{t) = 



^»+i,»+i(0 






(3.1) 



i = — 1. Here R{t) =|| r, •_,(<) ||. The 

operator /\* R{t) naturally acts on the i— th 
exterior power /\' i?” by the following rule: 
A‘ ^(0(^1 A ... A V,) = R{t)vi A ... R{t)vi for 
any Vi, • • • , u,- 6 We have, further, the fol- 
lowing relations: 



^ii(0 =< R{t)ei,R{t)ei >= 



< exp{0{0)t)ei^exp{0{0)t)€i > = 



Let 0{t) =11 Oij(t) II be the solution to the 
Toda flow 

0 = [0,7t0], 

such that 0(0) is upper Hessenberg orthogonal 
and irreducible. Here ttO = 0_ — Ol and 0_ 
is strictly lower triangular part of O. Then 
0{t) possesses the same properties and 0{t) 
converges when / — ^ oo to a block diagonal 
matrix. Each two by two block corresponds to 
a pair of complex conjugate eigenvalues. The 
blocks are arranged in the decreasing order of 
real parts of eigenvalues [7, 5]. From the pre- 
vious discussion we know that 0{i) is almost 



< ^i?L(^)ei >— Ai(<). 

And more generally 

I 

< ei A...Ae,-,/\r(<)(ei A...Ae,) >= A,(f), 

(3.2) 

i = 1, • • • , n. By (3.2) we easily obtain 

^«-H,;+l(0 _ V Aj+i(<)Aj-i(t) 

r.-,,(i) “ Aj(t) 

The result now follows by (3.1). ■ 
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We have the following differential equations for 

Gi : 

t = l,***,n — 1. Recalling that o,,,- = 
i = i/n = ±1, and i/f + af = 1, we 

obtain + a, a,- = 0 or 

CTi = 

i'i = - I't+i), 

i = — 1. It is interesting to find how 

moments r,(t) =< ei,0(t)‘ei > evolve under 
the Toda flow. Consider the family of rational 
functions 

ft{z) =< ei,[zE - 0(t)]"'e, > . 
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